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ABSTRACT 

Using 7- year data from the Wilkinson Microwave Anisotropy Probe I identify a sharp "edge" in the 
microwave haze at high Galactic latitude (35° < |6| < 55°) that is spatially coincident with the edge of 
the "Fermi Haze/Bubbles". This finding proves conclusively that the edge in the gamma-rays is real 
(and not a processing artifact), demonstrates explicitly that the microwave haze and the gamma-ray 
bubbles are indeed the same structure observed at multiple wavelengths, and strongly supports the 
interpretation of the microwave haze as a separate component of Galactic synchrotron (likely generated 
by a transient event) as opposed to a simple variation of the spectral index of disk synchrotron. In 
addition, combining these data sets allows for the first determination of the magnetic field within 
a radio bubble using microwaves and gamma-rays by taking advantage of the fact that the inverse 
Compton gamma-rays are primarily generated by scattering of CMB photons at these latitudes, thus 
minimizing uncertainty in the target radiation field. Assuming uniform volume emissivity, I find that 
the magnetic field within our Galactic microwave /gamma- ray bubbles is ~ 5 fiG above 6 kpc off of 
the Galactic plane. 

Subject headings: Galaxy: center — ISM: structure — ISM: bubbles — Radio continuum: ISM 



1. INTRODUCTION 

Recent full sky data sets by the Wilkinson Microwave 
Anisotropy Probe (WMAP) and the Fermi Gamma- Ray 
Space Telescope have revealed the presence of a new and 
very large structure in the Milky Way. This emission 
manifests as an excess of both microwaves and gamma- 
rays when removing Galactic diffuse emission associ- 
ated with known emission mechanisms at these wave- 
lengths and has come to be named the "WMA P Haze" 
dFinkbeinerl l200i iDobler fe Finkbeinerl [20081: IDoblerl 
1 20 12b and the " Fermi Haze /Bubbles" (jDobler et al .1120101 : 
ISu et al.l l2010f ) in the microwaves and gamma-rays re- 
spectively. 

In microwaves, the haze is synchrotron emission with 
a brightness temperature as a functio n of frequency 
T oc v$ H with p H w -2.5 (jDoblerl 120121 ) . This spectral 
dependence implies that the underlying electron spec- 
trum (number density as a function of energy) is given 
by dN/dE oc with 7^-2 for energies -10 GeV. In 
gamma-rays, the haze/bubbles is most likely due to in- 
verse Compton (IC) scattering of the starlight, infrared, 
and cosmic microwave background (CMB) interstellar 
radiation field (ISRF)0 In the original d i scove ry paper 
of the Fermi haze/bubbles, Dobl eret al.l (j201Qh showed 
that the spectrum of the gamma-rays is consistent with 
IC emission from the same electron population responsi- 
ble for the microwave haze synchrotron in both energy 
dependence a s well as ove rall normalization. However, 
as shown by IDoblerl ([2012), the existence of —1-10 GeV 
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3 There is also the possibility that the gamma-ray emission is due 
to the decay of 7r° particles generated by proton-proton collisions 
within the haze/bubbles. However, to match the amplitude and 
brightness profile, this sc enario relies on a ~ 10 9 yr wind (see 
Crocker & Aharonian 2011) and it is difficult to reconcile the sharp 
edge with this long timescale. In addition, there is no associated 
Ha signal as is typically seen in winds and the spectrum appears 
consistent with an inverse Compton scenario (see 331 . 



IC gamma-rays at high latitude where the ISRF is domi- 
nated by the CMB require electrons with energies ~TeV. 

Taken together, the microwaves and gammas imply 
that the electron spectrum is roughly a powerlaw from 
~ 1-1000 GeV suggesting that either the electrons have 
not had sufficient time to cool (the cooling time in this 
energy range is ~ 10 6 -10 7 yr; see ISu et al.ll2010l ) or are 
conti nuously being accelerated within the haze/b ubbles 
(e.g., IDobler "^FaI]l2QTll fMertsch fe Sa7^l20lTI) . The 
large volume of hard spectrum cosmic-rays has made 
it difficult to identify an underlying origin for the 
haze/bubbles and there have been num erous studies ex- 
plorin g the possibilities fro m starbursts (|Biermann et al.l 
2010), to Galactic winds (Crocker & Aharonian 2011), 
to jet blown bubbles (|Guo fc M athews 20lH; IGuo et all 
l201lh . and even co-a nnihilation of dark matter particles 
in the Galactic halo (jDobler et al.ll201lD . 

The goal of this letter is not to delve into the ori- 
gin question any further, but rather to address one 
of the main discrepancies between the microwave and 
gamma-ray haze / bubbles; namely , the observation that 
(as pointed out bv lSu et al.ll2010f ). the gamma-ray emis- 
sion appears to have a sharp edge at latitudes \b\ « 
50°, while th e microwaves fall off in intensity closer to 
\b\ « 35° (see IDoblerl l2012f ). In g2]I describe the compo- 
nent separation methods used to uncover both the mi- 
crowave and gamma-ray haze/bubbles in WMAP and 
Fermi data, and in $3] I will show that, in fact, the mi- 
crowaves also have a sharp edge at b ~ —50° that is 
spatially coincident with the edge in the Fermi emission. 
In 2] I summarize this work and describe the important 
implications for this microwave edge on the interpreta- 
tion of the Galactic haze/bubbles. 

2. METHODS 

The most straightforward method for uncovering the 
Galactic haze/bubbles in both the WMAP and Fermi 
maps is via "template fitting" which uses maps at other 
wavelengths to morphologically trace known emission 
mechanisms in the data. A simple linear regression of 
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Fig. 1. — The difference between the inferred and true spec- 
tral index as a function of the true index for the CMB estimate 
used in this letter. The inferred haze spect rum measured by 
IDobler & Finkbeinerl (120081 ) and IDoblerl (1201^ is (3 H « -2.55. 
Since the CMB bias for a T oc v$ foreground with (3 w —2.5 is 
only 0.015, the measured spectrum is not significantly biased. This 
implies that the electron population has an energy distribution of 
dN/dE oc E' 2 - 1 at E ~ 10 GeV. 

these templates (using an appropriate mask) against the 
data yields template amplitudes which can be used to 
"peel away" these foregrounds. In this letter, the tem- 
plates are fit at each frequency (energy) for the WMAP 
(Fermi) data implying that no constraint is put on 
the shape of the spectrum of each emission mechanism, 
though it is assumed that the spectrum does not vary 
significantly with position. The d etails of the template 
fitting used here can be found in IDobler fc Finkbeinerl 
(|2QQ8[ j and IDoblerl ([2012ft (hereafte r DFQ8 and D12 re - 
spectively) for the microwaves and IDobler et al.l (|2010ft 
for the gammas. 

2.1. WMAP analysis 

F or the WMAP analy sis, the templates used are 
the iSchlegel et all (Il998ft fSFD) dust map evaluated 
at 94 GHz by Finkbeiner_el_aL| (fl999h (FDS), the 
Hasl am et all (Jl982ft 408 MHz map, and the Finkbeiner 
(2003|) Ha composite map. These templates are meant 
to trace the three primary Galactic emission mechanisms 
at microwave wavelengths: thermal and spinning dust 
(FDS), soft synchrotron (Haslam), and free-free (Ha). 
In addition, I include the haze and hard disk bivariate 
Gaussian templates used in D12. Because the CMB is 
of comparable brightness to the haze at WMAP wave- 
lengths, I presubtract the CMB5 estimate for the haze 
given by DF08. Pixels for which the dust extinction at 
Ha is greater than 1 magnitude and for which the Ha 
intensity is greater than 10 Rayleigh are masked in the 
fit, as well all point sources in the WMAP and Planck 
ERCSC (30 GHz to 143 GHz) catalogs. 

These Galactic templates are fit to the WMAP data 
via the regression equation w u — c = Va v where w u 
is a map of the WMAP data at frequency c is 



the CMB5 estimate, and P is a matrix of template 
maps. The equation is solved for the coefficients a v = 
(P T n- 1 P)- 1 (P T n- 1 )(w i/ -c) where n is the mean noise 
in a given band. The haze residual in each band is de- 
fined as Rh = Wj/ — c — Va v + ajj x H where H is the 
haze template. That is, it is the residual of the regres- 
sion plus the amount of haze template removed. In order 
to account for spectral variations with position, the fit is 
performed independently on the regions given in D120 
A composite Rh is constructed from the union of those 
regions. 

2.2. Fermi analysis 

At Fermi energies, there are three main sources 
of diffuse gamma-ray emission from the Galaxy: 
bremsstrahlung, IC, and 7r° decay. The latter component 
is the dominant emission mechanism at high latitudes. 
For the purpose of comparing the haze/bubbles in Fermi 
to the emission in WMAP, this letter concentrates on 
the extreme high Galactic southern latitudes, b < —35°. 
High northern Galactic latitudes are contaminated by 
dust-correlated e mission in both the microwaves (likely 
spinn ing dust; see lDraine fc Lazari an 1998: IDobler et all 
2009) and the gammas (tt° decay) and make a clear sep- 
aration of the haze/bubbles emission difficult in both 
datasets. 

Because the 7r 0, s are created via collisions of cosmic- 
ray protons with the ISM, I use the SFD map of dust 
column density as a template for this emission since it is 
a reasonable tracer of dust and gas in our Galaxy. It is 
important to note that this template is an integrated col- 
umn density while 7r° emission is proportional to the ISM 
density times the proton number density, and so line-of- 
sight effects render the SFD map an imperfect tracer of 
7T° emission. This has dramatic effects at lo w latitudes in 
the ha ze /bubbles regions as pointed out bv IDobler et al.l 
(|2011l ). however for b < —35° as is used here , a simple two 
temp late (SFD plus bivariate Gaussian; see IDobler et al.l 
l2Q10h description of the data is sufficient to isolate the 
Fermi h aze /bubbles emissio n. 

As in IDobler et al.l (|2010ft the fit minimizes the log- 
likelihood, ln£ = ^Jfciln/ii —lii — ln(/^!)], where ki is 
the map of observed counts at pixel z, fi is a synthetic 
counts map given by ja(E) = S(E) x (mask) x (exposure), 
and S(E) is the synthetic sky map S(E) = 6sfd x SFD + 
fruni- This log-likelihood is minimized over the 
template amplitudes &sfd, and b un { (a spatially uni- 
form contribution) for pixels outside of the mask and for 
b < — 35 °0 I emphasize that this is not meant to rep- 
resent a perfect model for the diffuse emission from our 
Galaxy, but rather serves to effectively isolate the Fermi 
haze/bubbles emission so that it can be morphologically 
compared to the microwave emission. For a thorough 
anal ysis of Galactic diffuse emission obse rved by Fermi 
see iThe Fermi-LAT Collaboration! (j2012ft . Lastly, the 
analysis was also redone using the Fermi diffuse model 
(gll_iem_v02 . f it0) as opposed to the SFD map as a 

4 An extreme example of this technique are pixel-by-pixel fits of 
the data which use a combination of spectral and spatial templates. 
While the flexi bility of th e se m odels makes the haze analysis more 
difficult, Pietrobon et al. (2011) showed that the microwave haze 
is indeed recoverable with these techniques. 

5 All gamma-ray results in this letter are derived with 1. 6-year 
Fermi maps constructed as described in IDobler et al. ( 2010). 

^ http : //f ermi . gsf c . nasa . gov/ssc/dat a/access/lat /BackgroundModels . html 
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Fig. 2. — Top row: the microwave haze at WMAP K- and Ka-band smoothed to 1°. As noted in Dobler (201.2), the microwaves seem to 
fade quickly below b ~ —35° in contrast to the Fermi haze/bubbles (lower left) which continue down to b « —50°. The Fermi residual is 
just the difference between the data and the Fermi diffuse model for visualization (see However, when smoothing the WMAP data to 
2° and stacking with weights given by (V/^k) 2 ' 5 ^ can ^ e seen m the lower right panel that there is an edge in the microwaves at b ~ —50° 
as well. 



template and the results are not significantly changed. 
3. RESULTS 

Directly comparing the haze/bubbles residuals at high 
latitude will involve weighted stacking of Rh in multi- 
ple WMAP bands, and so it is useful to estimate the 
extent to which the spectrum of the microwave emission 
is affected by the CMB bias described in DF08. This 
bias comes from the fact that any CMB estimate used in 
$2] will inherently have some residual foregrounds after 
cleaning emission from the Galaxy, and since that CMB 
estimate is presubtracted with a fixed CMB spectrum, 
this imprints a bias on the inferred haze spectrum (see 
DF08). The CMB estimate used here is generated by 
the linear combination of the WMAP data that mini- 
mizes the variance in unmasked pixels of c = . Cj w j 
where w^- is the WMAP data minus the FDS prediction 
for thermal dust emission in each band j and . Q = 1 
for w in thermodynamic AT units. 

This last constraint on £ has the consequence that the 
measured spectrum cannot be used to infer the true spec- 
trum. That is, while the amplitude of the bias can be 
estimated as in DF08, the exact bias cannot be known. 
However, it is possible to answer the question: for a 
foreground with a given true spectrum T anten na oc iA rue , 
what would be the inferred spectral index /3i n f er ? This 

5j3 = ftnfer — Ptrue I s shown as a function of Ptrue m 

Figure HJ Despite the very large possible CMB bias, 
for a true spectral index of f3 true = —2.5, the inferred 
spectral index when comparing K- to Ka-band would 
only be biased by 5f3 = —0.015 given the CMB5 coeffi- 



cients C = (0.11, -0.48, 0.12, 0.16, 1.09). The implication 
is that the spectrum measured by DF08 and D12 for the 
microwave haze Tu oc v$ H with fin ~ —2.5 is likely not 
significantly biased and the underlying electron popula- 
tion is very close to dN/dE oc E~ 2 . 

The full sky residuals and — the haze/bubbles 
in WMAP and Fermi data — are shown in Figure [2] 
and illustrate several of the key characteristics noted 
in previous studies: scaling Rjj by v 2 ^ yields roughly 
equal brightness ind i cating a T oc v~ 2 - 5 spectrum 
(jDobler fc Finkbeinerl l200l ). B% a nd are roughly 
spatially coincident at low latitudes (jDobler et al.l 1201(1 
ISu et aLl l2Q10h , at 1° s moothing R% appears to fade 
quickly below b = -35° (jDobl er 2012), and i ?£ appears 
to have a sharp "edge" at \b\ ~ 50° (jSu et aLll2010[ ). The 
lack of a similar edge in the WMAP data at b ~ —50° 
in WMAP has led to some ambiguity about whether the 
WMAP haze and the Fermi haze/bubbles are in fact the 
same structure. However, creating the weighted stack 
of the WMAP emission R% + (^Ka/^x) 2 - 5 x and 
smoothing to 2° (i.e., the same smoothing as the Fermi 
map) reveals an edge in the microwave haze that appears 
coincident with the southern edge in the gamma-ray bub- 
bles. 

To assess the significance of this feature, I zoom in on 
the extreme southern latitudes with b < —35° in Figure 
[3] (corresponding to heights > 6 kpc above the GC), and 
a clear sharp edge is evident in b oth the gamma -rays 
and the microwaves. As noted in ISu et all ([2010) the 
center of the Fermi bubble at these latitudes is roughly 
i rsj —4.5°. Binning the sky into polar bins centered 
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Fig. 3. — Top left: the southern Fermi bubble at latitudes b < —35° where the edges are most discernible in the gammas. Top right: 
the same region but for the WMAP K-band haze/bubbles smoothed to 2°. There is a clear spatial correspondence between the Fermi 
haze/bubbles and the WMAP haze, including an edge at high latitudes. Taking the center of the emission to be (£, b) C en = (—4.5°, —35.0°) 
(dashed white line) and plotting the intensity of the gammas and microwaves as a function of distance from the center (i.e., integrating 
over the angular bins shown in dotted lines in the upper right panel) reveals an unambiguous detection of an edge at high latitudes that is 
spatially coincident with the edge in the Fermi haze/bubbles as shown in the bottom left panel. This emission is well fit by a flat brightness 
profile with a sharp edge at r ~ 17° from the center (performing the same annular bin using £ cen = 180° yields no such feature). Lower 
right: the same but integrating over radial bins (for two ranges in r) as a function of angular bin. 

the bubble center. For the 2° smoothing shown here, 
the statistical significance of the edge identified in this 
way is high for both the K-band microwave haze/bubbles 
and the K-, Ka-band weighted stack (a null test calcu- 
lated by setting (4en, ^cen) = (—180°,— 35°) shows no 
evidence for an edge at the Galactic anti-center). Finally, 
by integrating a uniform brightness, infinitely sharp edge 
(smoothed to 2° and with the same mask) in the same 



on the bubble center (4en,^cen) = (—4.5°, —35°), inte- 
grating annuli with unmasked pixel latitudes < —35° 
(see figure), and plotting as a function of distance from 
the bubble center, the lower left panel of Figure [3] shows 
the so uthern Fermi bubble edge described in ISu et al.l 
(2010). The same plot generated with the microwaves 
also shows a clear edge that is spatially coincident with 
the Fermi bubble edge, located at roughly r = 17° from 



5 



CO 
CO 

CM 

E 

o 
> 



1.8 
1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 



23 GHz microwaves 

23-33 GHz stacked 

l . _ flat profile 

T T[r 23GHz,l cen = 180° 



2-5 GeV gammas 



upper curves, r < 20° 
lower curves, r > 20° 




32.6 



18.91 



10 20 30 
Radial bins [deg] 



50 100 150 

Angular bins [clockwise, b<-35°] 



Fig. 4. — The same as the bottom panels of Figure|3] but using the Fermi diffuse model in the template fit. 
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Fig. 5. — The spectrum of the Fermi haze/bubble for b < —35° 
and within r < 17° of (£,b) = (-4.5°, -35°). The four different 
lines represent spectra derived from haze template or residual am- 
plitudes for both SFD and Fermi diffuse model fits. The spectra 
have been normalized to the SFD residual spectrum at 4 GeV. 
While the uncertainty at low energies is large, the four spectral 
estimates agree above 4 GeV and are consistent with IC emission 
from an electron population with dN/dE oc E~ 2 . The spectrum is 
softer than at lower latitudes (suggesting that the IC signal is due 
primarily to CMB scattering at distances > 6 kpc above the Galac- 
tic plane), which argues for a leptonic (IC) rather than hadronic 
(7T° decay) origin. 

way, Figure [3] shows that Fermi and WM AP bubbles are 
consistent with a sharp edge. 

While this is a clear detection of an edge in the WMAP 
haze, there is angular dependence in the microwave emis- 
sion that is not readily apparent in the gamma-rays. Fig- 
ure 02 shows the emission as a function of annulus angle 
both inside and outside the haze/bubbles (r < 20° and 
r > 20° from (4 en , b cen )) in the lower right panel and in- 



dicates an "arm" of emission in the microwaves for annu- 
lar angles less than ~ 30°. No such arm is evident in the 
gamma-rays. For annular angles greater than 30° there is 
no other clear structure present in the microwaves though 
the microwave bubble emission is detected in each an- 
gular bin (i.e., there is an excess of emission interior 
to the microwave edge compared to exterior). Figure 
H] shows the same results, but using the Fermi diffuse 
model template in the fit. While the overall amplitude 
of the gamma-ray signal is slightly lower, the edge is still 
spatially coincidence with the edge in microwaves. 

Figure [5] shows the spectrum of the Fermi haze/bubble 
for b < —35° and within r < 17° of the bubble center 
for four different spectral estimates: using the SFD map 
or the Fermi diffuse model in the template fit and esti- 
mating the spectrum via unmasked pixels in x H 
or Rh- Below ~ 4 GeV, the four different types of 
spectral estimates differ, indicating a systematic bias in 
the derived spectrum due to imperfect template approxi- 
mations and/or brighter, softer components leaking into 
the residuals. However, above ~ 4 GeV, the spectra 
are all consistent with IC from electrons scattering pre- 
dominantly off of CMB photons and having a spectrum 
dN/dE oc E~ 2 (i.e., the same as that required to make 
the T oc v~ 2 ' b microwave haze spectrum). It is interest- 
ing to note that this spectrum of the "cap " of the bub- 
ble is some what softer tha n that found by iDobler et al.l 
(|2010D and ISu et all (|2010l ) who include lower latitudes, 
in agreement with the IC scenario (since the higher en- 
ergy optical and IR components of the ISRF have lower 
amplitude at high latitudes j3 and inconsistent with a 
hadronic scenario, which is also disfavored by the de- 
tection of the microwave haze at b < —35° in Figure [H 

Since the IC emission is a function of the number den- 
sity of electrons dN/dE and the ISRF intensity while 
the synchrotron intensity depends on dN/dE and the 
magnetic field an estimate of B can be made un- 
der several simplifying assumptions. Given the (CMB 

7 For this calculation, the GALPR0P ISRF estimate was used (see 
http://galprop.stanford.edu/). Since this estimate only goes 
up to ~ 5 kpc in height above the plane, I take the optical and IR 
components to be 75% of the 5 kpc value. 
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dominated) ISRF model above and assuming the same 
dN/dE = No x E~ 2 for both signals, a uniform mag- 
netic field and volume emissivity model above 6 kpc 
(b < —35°) yields B ~ 5 fiG within the southern bubble. 



4. SUMMARY 

Using 7-year WMAP data to isolate the microwave 
haze and comparing this to the Fermi haze/bubbles at 
southern Galactic latitudes less than —35°, I have pre- 
sented the detection of a sharp edge in the microwave 
haze that is coincident with the Fermi bubble edge. This 
microwave bubble edge is evident when smoothing the 
haze to 2° in both the K-band WMAP data as well 
as a stack of K- and Ka-band weighted by ^ 2 5 , where 
v is the microwave band. I have also shown explicitly 
th at, for the CM B estimate used in this study as well 
as iDoblerl (j2012l ). the u~ 2 - b spectrum of the microwave 
haze/bubbles is not significantly biased by systematics 
indicating that the electrons responsible for generating 
this synchrotron have a number density as a function of 
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energy dN/dE oc E~ 2 , in excellent agreement with the 
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The detection of an edge in the microwave haze ~ 
50° above the plane and coincident with the Fermi 
haze/bubbles has several important consequences. First, 
it proves conclusively that the microwave and and the 
Fermi bubbles are the same structure observed at mul- 
tiple wavelengths. Second, given the vastly different ex- 
periments, the detection of an edge in microwaves proves 
that the edge in gammas is real and not due to process- 
in g artifacts or low pho ton counts in the m aps generated 
by iDobler etall (j2Q10D and lSu et all (|2Q10f ). Finally, the 
sharp edge, coupled with the hard spectrum of the emis- 
sion, suggests a transient event for the origin of the mi- 
crowave haze indicating that it is a separate component 
of diffuse emission in our Galaxy and not merely a spatial 
variation in the spectral index of the disk synchrotron. 
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